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We investigate the asymptotic relaxation of so-called pulled fronts propagating into an unstable 
state. The "leading edge representation" of the equation of motion reveals the universal nature of 
their propagation mechanism and allows us to generalize the universal algebraic velocity relaxation 
of uniformly translating fronts to fronts, that generate periodic or even chaotic states. Such fronts 
in addition exhibit a universal algebraic phase relaxation. We numerically verify our analytical 
predictions for the Swift-Hohenberg and the Complex Ginzburg Landau equation. 
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Many systems, when driven sufficiently far from equi- 
librium, spontaneously organize themselves in coherent 
or incoherent patterns. The ubiquity of such structures in 
almost all fields of the natural sciences has inspired 
much of the recent scientific effort to uncover the var- 
ious mechanisms underlying their behavior. Especially 
in physics, the insight from the seventies that in crit- 
ical phenomena universality classes are determined es- 
sentially by the symmetry of the order parameter and 
the dimensionality, initially raised some hopes that there 
would be analogous broad universality classes in noncqui- 
librium pattern formation. Over the last two decades, it 
has become clear, however, that such far-reaching uni- 
versality does not exist: While there are various general 
dynamical and instability mechanisms, there is not al- 
ways a sharp selection mechanism. Moreover, if there is 
sharp selection, the particular mechanism may depend on 
the specific boundary conditions, initial conditions, etc. 

In the case of front propagation, there have, neverthe- 
less, been several hints of a generic dynamical mecha- 
nism There is a large class of fronts propagating 
into an unstable state whose asymptotic velocity equals 
V* , the asymptotic spreading velocity of linear pertur- 
bations about the unstable state. Such fronts are called 
pulled fronts, as they are "pulled along" by the leading 
edge of the profile whose dynamics is governed by the 
linearization about the unstable state It is the 
purpose of this letter to show that within the subclass 
of pulled front propagation, a remarkable degree of uni- 
versality does hold: Irrespective of whether such fronts 
are uniformly translating or generate periodic or chaotic 
patterns, the velocity v{t) and phase T{t) of pulled fronts 
which emerge from steep initial conditions (falling off 
faster than e~^ ^ for x oo), display a universal power 
law relaxation with time t, expressed by 

v{t) = v*+X{t) (1) 



^ (3) 

As explained below, the coefficients v* ,k* — q*+iX* , and 
D are all given explicitly in terms of the dispersion rela- 
tion of the linearized equation. We shall focus on deter- 
mining how these exact asymptotic relaxation formulas 
emerge, and why they are independent of the nonlineari- 
ties, the precise initial conditions, or on whether the front 
dynamics is regular or chaotic. Before embarking on this, 
however, it is important to explain what we mean by ve- 
locity and phase for the various types of fronts. 

Uniformly translating pulled fronts. The simplest types 
of fronts are those for which the dynamical field 4>{x, t) 
asymptotically approaches a uniformly translating pro- 
file (j) = ^v'iO: i = x — v*t, as happens, e.g., in the cele- 
brated nonlinear diffusion equation dt(j) = d'^4'+4'~4'^ for 
fronts propagating into the unstable = state. If we 
define level curves as the lines in an x,t diagram where 
(j){x, t) has a particular value, we can define the velocity 
v{t) as the slope of a level curve. For uniformly trans- 
lating fronts, q* —0 = ImZ3; (|^) then reduces to the ex- 
pression derived for uniformly translating fronts in [Q. 
The remarkable point is that the expression for v{t) is 
in this case completely independent of which level curve 
one traces. Moreover, it was shown in Q that the nonlin- 
ear front region is slaved to the leading edge of the front 
whose velocity relaxes according to (^). This results in 

(t>{x,t)^<^,i^^){ix) + 0{t-^) , ^x^Vt, (4) 
U = x-v*t-X{t) , (5) 

where ^ — x — vt solves the o.d.e. for a front prop- 

agating uniformly with velocity v. v{t) in (^ is the in- 
stantaneous velocity of the front, and the frame is 
shifted by the time dependent quantity X{t). Since the 
collective coordinate X{t) diverges as \nt for large t ac- 
cording to (^ , the difference between ^x and a uniformly 
translating frame is crucial — only in the former can we 
follow the relaxation. Uniformly translating fronts have 
no phase, hence all terms in (0) vanish identically. 
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Coherent pattern generating fronts. As an example of 
coherent pattern generating fronts, we consider the so- 
called Swift-Hohenberg (SH) equation 



dtu = £w - (1 + dD'^u 



£ > 



(6) 



The space-time plot of Fig. 1(a) illustrates how SH- 
fronts with steep initial conditions (falling off faster than 
^ as x^oo into the unstable state u = Q) generate a 
periodic pattern. It is known that they are pulled 
In this case, new level curves in an x, t plot are constantly 
being generated. If we define in this case the velocity as 
the slope of the uppermost level curve, one gets an os- 
cillatory function. Its average is v{t) given in (|l|), but 
v{t) is difficult to extract this way. Numerically, it is 
better to determine the velocity from an empirical enve- 
lope obtained by interpolating the positions of the max- 
ima. Since these pattern forming front solutions for long 
times have a temporal periodicity u(^, t) = m(^, t-\-T) in 
the frame = x — vt moving with the velocity v of the 
front, the asymptotic profiles can be written in the form 
E„=i e-2'r"*/^f7;;'(C) + c.c. In terms of these complex 
modes U, our result for the shape relaxation of the pulled 
front profile becomes in analogy to (H) 



u{x,t) ~ ^( 



(7) 



n=l 



with the frequency fi* given below [|To[. Eq. (0) shows 
that T(t) is the global phase of the relaxing profile, as the 
functions only have a ^jf-dependence. The result of 
our calculation of the long time relaxation of v{t) and 
r(i) is given in (^) - (^). We stress that while for e ^ 0, 
an ansatz like (p ) leads to an amplitude equation for the 
n — ±1 terms, our analysis applies for any e > 0. 

Incoherent or chaotic fronts. The third class we con- 
sider consists of fronts which leave behind chaotic states. 
They occur in some regions of parameter space in the 
cubic Complex Ginzburg-Landau equation ||ll|] or in the 
quintic extension (QCGL) that we consider here, 

dtA = + (1 + iCi)dlA + (1 + iC^)\A\''A 

-{l-iC^)\A\'A. (8) 

Fig. 1(h) shows an example of a pulled front in this equa- 
tion. Level curves in a space-time diagram can now also 
both start and end. If we calculate the velocity from the 
slope of the uppermost level line, then its average value is 
again given by (||) but the oscillations can be quite 
large. However, our analysis confirms what is already 
visible in Fig. 1(b), namely that even a chaotic pulled 
front becomes more coherent the further one looks into 
the leading edge of the profile. Indeed we will see that in 
the leading edge where |A| ^ 1 the profile is given by an 
expression reminiscent of (|^), 

A{x, t) « e-*"**-*i^We*'=*«^ V(ex), 1 « U « Vt. (9) 



The fluctuations about this expression become smaller 
the larger ^x- 

In Figs. 1(c) and 2(c) we show as an example results 
of our simulations of the SH equation (|q) and the QCGL 
(^ . They fully confirm our predictions ^) and (||) for the 
asymptotic average velocity and phase relaxation. Note 
that for the QCGL, the fluctuations are indeed smaller 
the more one probes the leading edge region. 

We now summarize how these results arise. 

Calculation of the asymptotic parameters. Although 
this is well-known we first briefly summarize how 
the linear spreading velocity v* and the associated pa- 
rameters A* etc. arise, as the analysis also motivates the 
subsequent steps. After linearization about the unsta- 
ble state, the equations we consider can all be written 
in the form dt(f) — £(9a;,9^, • • ■)(j). For a Fourier mode 
^-lujt+tkx ^ this yields the dispersion relation ui{k). The 
linear spreading velocity v* of steep initial conditions is 
then obtained by a saddle point analysis of the Green's 
function G of these equations. In the asymptotic frame 
^ = a; — v*t, G(^, t) becomes 



„ik'i-iQ.*t ^ 



4Dt 



y/AnDt 



for large times. Here 51(fc) = uj{k) — v*k, and 



dflik) 



dk 



= , lmn{k*) = , 



id^n{k) 



2d/fc2 



(10) 



(11) 



The first equation in (|11|) is the saddle point condition, 
while the second one expresses the self-consistency con- 
dition that there is no growth in the co-moving frame. 
These equations straight forwardly determine v*,k* ~ 
q* + iX* , D and the real frequency ft* — Q{k*) |]l5|] . 

Choosing the proper frame and transformation. Eq. 
( |lo| ) not only confirms that a localized initial condition 
will grow out and spread asymptotically with the velocity 
V* given by (|ll]), but the Gaussian factor also determines 
how the asymptotic velocity is approached in the fully lin- 
ear case. Our aim now is to understand the convergence 
of a pulled front due to the interplay of the linear spread- 
ing and the nonlinearities. The Green's function expres- 
sion ( |lo|) gives three important hints in this regard: First 
of aU, G{^,t) is asymptotically of the form 6'*^*?^'^'* 
times a crossover function whose diffusive behavior is 
betrayed by the Gaussian form in (Rw. Hence if we 
write our dynamical fields as A — e^^ 'V'(? 1 1) for the 
QCGL (|) or u = e''=*«-^*^**V(C, t)+c.c. for the real field u 
in (|^), we expect that the dynamical equation for V'(Ci 
obeys a diffusion-type equation. Second, as we have 
argued in [Q, for the relaxation analysis one wants to 
work in a frame where the crossover function ■0 becomes 
asymptotically time independent. This is obviously not 
true in the ^ frame, due to the factor l/\/t in (px|). 
However, this term can be absorbed in the exponential 
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prefactor, by writing ^-^ g^^'f-^f^** = ^ik'^-in't-uint ^ 
Hence, we introduce the logarithmically shifted frame 
S,x ~^—X{t) as already used in (||). Third, we find a 
new feature specific for pattern forming fronts: the com- 
plex parameters, and D in particular, lead us to introduce 
the global phase r{t). We expand t{t) like X{t) § 



Xit) 



ci C3/2 

t ^ ^3/2 + 



t 



"3/2 



(12) 



and analyze the long time dynamics by performing a 
"leading edge transformation" to the field i/', 

QCGL: A = f.^k■'(,x-^irt-^T(t) ^^^^t) , 

SH: U = f^k'ix-^^Vt-^T(t) ^(^^^ + c.c. (13) 

Steep initial conditions imply that '(/'(Cx,^)— >0 as £,x^ 
00. The determination of the coefficients in the expan- 
sions (|l|) of X and f is the main goal of the subsequent 
analysis, as this then directly yields Eqs. (||) and (^). 

Understanding the intermediate asymptotics. Substi- 
tuting the leading edge transformation ( |l3|) into the non- 
linear dynamical equations, we get 

n— 3 

+ [i:(t)(% + ik*) + ir{t)]^p - iv(^) , (14) 

with Dn — (— i/n!)d"i7/(diA;)"|fc. the generalization of 
D in (|ll] ) (of course, for the QCGL, n{k) is quadratic 
in k, so Dn — 0)- In this equation, N accounts for the 
nonlinear terms; e.g., for the QCGL, we simply have 



N = e-^^'^^'lipl^ip [l~iC3 + {l~iC5)e 



The expression for the SH equation is similar. 

The structure of Eq. (|l^) is that of a diffusion-type 
equation with l/t and higher order corrections from the 
X and r terms, and with a nonhnearity N. The crucial 
point to recognize now is that for fronts, N is nonzero 
only in a region of finite width: For — > 00, N decays 
exponentially due to the explicit exponential factors in 
(|l5|). For £_x^—oo, N also decays exponentially, since u 
and A remain finite, so that "0 decays as accord- 
ing to (|l3|). Intuitively, therefore, we can think of (14) as 
a diffusion equation in the presence of a sink N localized 
at some finite value of ^x- The ensuing dynamics of tp 
to the right of the sink can be understood with the aid 
of Figs. 2(a) and (b), which are obtained directly from 
the time-dependent numerical simulations of the QCGL 
(^. To extract the intermediate asymptotic behavior il- 
lustrated by these plots, we integrate ( p^ ) once to get 

dt f^^d^x yj^Dd^^^p + Y, d^^'^p + (16) 

+i[k*xit)+t{t)] f^^dex^+x{t)^- f^^dexNiij) 



Now, in the region labeled I in Fig. 2(b), we have for fixed 
^x and t ^ 00 that the terms proportional to X and F 
can be neglected upon averaging over the fast fiuctua- 
tions; the same holds for the term on the left. Since the 
integral converges quickly to the right due to the expo- 
nential factors in TV, we then get immediately, irrespec- 
tive of the presence of higher order spatial derivatives 



hm D—- 



d^xNii;) = aD 



(17) 



Here, the overbar denotes a time average (necessary for 
the case of a chaotic front). Thus, for large times in re- 
gion I, « a^x -l- /? in dominant order. Moreover, from 
the diffusive nature of the equation, our assertion that 
the fiuctuations of tp rapidly decrease to the right of the 
region where N is nonzero comes out naturally. In other 
words, provided that the time- averaged sink strength a 
is nonzero, a ^ 0, one will find a buildup of a linear gra- 
dient in \tp\ in region I, independent of the precise form 
of the nonlinearities or of whether or not the front dy- 
namics is coherent. This behavior is clearly visible in 
Fig. 2(b). We can understand the dynamics in regions 
II and HI along similar lines. In region HI the dominant 
terms in ( p^ are the one on the left and the first one on 
the second line, and the cross-over region II which sep- 
arates regions I and III moves to the right according to 
the diffusive law £_x ^ D\/t. 

Systematic expansion. These considerations are fully 
corroborated by our extension of the analysis of An- 
ticipating that "0 falls off for ^ 1, we split off a Gaus- 
sian factor by writing tp{S^x,t) — G{z,t) e~^ in terms of 
the similarity variable z = ^x'^ /(ADt), and expand 

G(z, t) = Vi (z) + 9o{z) + t-^/^gi (z) + • • • (18) 



This, together with the expansion (|l^) for X(t) and 
r(i), the left "boundary condition" that ip{(^x,t^oo) = 
(^£,x +/3 and the condition that the functions g{z) do not 
diverge exponentially, then results in the expressions (|^) 
for X{t) and (|) for f §. For the QCGL, the analysis 
immediately implies the result (|^) for the front profile in 
the leading edge. In addition for the SH equation, one 
arrives at (^ for the shape relaxation in the front interior 
along the lines of Q: Starting from the o.d.e.'s for the 
t/", one finds upon transforming to the frame that to 
0(t^^), the time dependence only enters parametrically 
through v{t). This then yields (^. 

In conclusion, we have shown that the long time relax- 
ation of pulled fronts is remarkably universal: Indepen- 
dent of whether fronts are uniformly translating, pattern 
generating or chaotic, the velocity and phase relaxation is 
governed by one simple formula, with universal dominant 
and subdominant power law expressions. 
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FIG. 1. (a) Space-time plot of a pulled front in the SH eq. ^ with e — 5 and Gaussian initial conditions. Time steps between 
successive lines are 0.1. (b) A pulled front in the QCGL eq. with e = 0.25, Ci = 1, C3 = Cs = —3, and Gaussian initial 
conditions. Plotted is |yl(a;,t)|. Time steps between lines are 1. (c) Scaling plot of the velocity relaxation {v{t) — v") ■ Tv/\ci\ 
vs. 1/r with r = t/T^ and characteristic time = [0^/2/^1)'^ ■ Plotted are from top to bottom the data for the SH eq. for 
heights u = 0.0001^, O-Ol^i, and (e = 5) as dashed lines, and for the QCGL eq. (|) for heights \A\ = 0.00002, 0.0002, and 
0.002 as dotted lines. The solid line is the universal asymptote —1/t + 
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FIG. 2. (a) and (b): Simulation of the QCGL eq. as in Fig. 1 (b) for times t = 35 to 144. (a) shows \N\ ( |l5| ) as a function 
of £,x- (b) shows I?/) I, which in region I builds up a linear slope 7/; cx Q^x, and in region III decays 1 ike a Gaussian widening 
in time. The lines in region II show the maxima of ip{^x,t) for fixed t and their projection ~ Vt into the {^x,t) plane. 
(c) shows the scaling plot for the phase relaxation. From top to bottom: SH (dashed) for u — 0.0001-^6, COl-^e, and ^/e 
(e = 5), and QCGL (dotted) for \A\ = 0.002, 0.0002, and 0.00002. Plotted is r{t) ■ Tr/ci vs. 1/r. Here r = t/Tr, and 
Tr = T„ ■ [1 + A*ImZ)"^''^/(q*ReZ)"^''^)] . The solid line again is the universal asymptote -1/r + l/r^^^. 
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